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TECHNICAL  MEMORANDUM  I598 


FORMATION  OF  A  VORTEX  AT  THE  EDGE 

OF  A  PLATE^ 

By  Leo  Anton 

1.  INTRODUCTION 


If  a  plate  is  moved  from  the  state  of  rest  transversely  to  its 
width  through  a  fluid  at  rest,  a  small  vortex  forms  at  every  edge  at 
the  start  of  the  motion;  'by  continual  added  infliix  of  new  parts  of  the 
fluid,  that  vortex  becomes  larger  and  grows  into  the  flow.   In  case  the 
flow  is  perpendicular  to  the  plate,  the  position  of  the  two  vortices  is 
symmetrical  at  the  beginning  of  the  motion  of  the  fluid.   In  a  later 
stage  of  development  it  is,  as  von  Karman  has  proved,  unstable.   Then 
new  vortices  originate  alternately  at  the  upper  and. lower  end  of  the 
plate  which  group  themselves  behind  the  plate  in  a  certain  manner  (von 
Karman 's  vortex  street). 

In  what  follows,  we  shall  treat  the  initial  flow  just  described  for 
a  plate  perpendicular  to  the  approaching  flow,  where,  therefore,  flow 
direction  and  plate  form  a  right  angle .   Since  in  this  case  the  vortices 
are  in  the  positions  of  reflected  images,  our  investigation  is  limited  to 
the  formation  and  the  growth  of  one  vortex  at  one  plate  edge.  We  visualize 
the  plate  as  suddenly  set  in  motion  and  then  moved  uniformly  at  constant 
velocity.   For  a  system  of  axes  moved  with  the  plate,  the  velocity  at 
infinity  is  therefore  to  remain  invariable  throughout  the  flow  duration. 
Moreover,  we  assume  the  plate  to  be  laterally  extended  to  infinite  length 
so  that  we  deal  with  a  plane  nonstationary  flow  phenomenon. 

Solution  of  the  problem  is  considerably  facilitated  by  disregard  of 
the  friction.   This  is  permissible  since  the  friction  of  most  fluids 
coming  into  question  is  very  small.   Thus  we  take  a^  a  basis  an  ideal 
fluid  where,  therefore,  the  viscosity  is  zero,  and  require  furthermore 
that  it  be  homogeneous  and  incompressible.   In  this  fluid  the  flow 
described  is  then  to  be  interpreted  as  a  potential  flow  with  free  vor- 
tices distributed  over  a  discontinuity  surface,  a  so-called  vortex  sheet. 
The  velocity  field  of  such  a  flow  is  uniquely  determined  by  the  velocity 
at  infinity  and  by  the  form  and  vortex  distribution  of  this  surface. 

A  solution  for  the  motion  of  a  fluid  about  a  plate,  starting  from 
the  state  of  rest,  has,  so  far,  been  achieved  only  for  infinitely  small 
angles  of  attack.   This  borderline  case  has  been  treated  by  H.  Wagner 

Ausbildung  eines  Wirbels  an  der  Kante  einer  Platte,"  Gottlnger 
Dissertation,  Ingenieur-Archiv,  vol.  X,  1939^  PP-  411-^4-27. 
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(ref.  1).   Here  the  free  vortices  lie  in  a  plane  vortex  sheet  which 
extends  from  the  trailing  edge.   The  solution  of  the  flow  problem  consists 
in  the  calculation  of  the  distribution  density  of  the  discontinuity  sur- 
face. Wagner  performed  these  calculations  for  accelerated  and  for  uniform 
motion  of  the  plate.  For  finite  angles  of  attack,  as  in  our  case,  the 
discontinuity  surface  forms  a  strongly  curved  vortex  sheet  which  obviously 
is  rolled  up  into  a  spiral  vortex.   Thus,  besides  the  distribution  density, 
the  shape  of  the  discontinuity  surface  also  is  unknown. 

We  now  make  the  following  simplification.   At  the  very  beginning  of 
the  motion  of  the  fluid,  the  magnitude  of  a  vortex  at  the  plate  edge  in 
proportion  to  the  width  of  the  plate  is  still  very  small.   Consequently 
the  vortices  at  both  edges  still  lie  outside  of  their  mutual  interference 
domain  so  that  one  deals  practically  with  a  flow  about  a  plate  of  semi- 
infinite  width.   This  means  a  great  advantage.   Since  no  unit  length 
exists  which  would  influence  the  process,  the  flow  conditions  must  be 
independent  of  the  magnitude  of  the  vortex  at  the  time .   The  growth  of  the 
vortex  therefore  consists  in  a  similar  increase  of  it.   In  this  connection 
the  treatises  of  L.  Prandtl  (ref.  2)  and  H.  Kaden  (ref.  5)  must  be  men- 
tioned.  Prandtl  performed  a  general  investigation  concerning  the  flow 
around  a  corner  for  certain  laws  of  acceleration.   No  solution  for  the 
plate  flow  resulted  from  the  Prandtl  formulation;  however,  one  can  derive 
from  it  the  laws  of  similitude.   In  the  treatise  mentioned,  Prandtl  sur- 
mises that  the  rolled-up  vortex  sheet  in  its  interior  may  be  of  the  type 
of  a  spiral  R  =  const. /cp"^   (R  and  cp  signify  polar  coordinates  counted 
from  the  center  of  the  spiral,  m  is  a  number).   Kaden  treats  the  rolling 
up  of  an  unstable  discontinuity  surface,  unilaterally  of  infinite  length, 
into  a  spiral  vortex.   As  we  shall  see  later,  Kaden 's  problem  becomes 
identical  with  ours  for  the  interior  of  the  vortex.   The  same  laws  of 
similitude  are  valid  for  both  cases.   The  solution  for  the  vortex  core 

/  2/5 
R  =  const. /cp  '    found  by  Kaden  applies,  therefore,  also  to  the  interior 

of  our  vortex  sheet;  the  above-mentioned  prediction  of  Prandtl  is  thereby 

confirmed. 

The  outer  part  of  the  discontinuity  surface,  thus  the  transitional 
region  from  the  spiral-shaped  vortex  core  to  the  plate  edge,  we  procure 
with  the  aid  of  graphical  methods.   The  solution  obtained  is  then  used 
for  finding  the  flow  about  the  plate  of  finite  width. 


2.  OUTLINE  OF  THE  METHOD 


In  this  section  we  shall  briefly  describe  the  method  by  which  one 
arrives  at  solution  of  the  problem.  We  shall  start  from  the  plate  of 
finite  width,  and  shall  then  treat  the  plate  of  infinite  width  as  a 
special  case. 
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In  the  case  of  the  vortex-free  potential  flow  about  the  plate  placed 
transversely,  the  velocity  at  the  edges  is  known  to  become  infinitely 
large;  however,  under  actual  conditions  this  is  never  possible.  We  must 
therefore  demand  also  in  this  case  that  the  velocity  there  always  has 
finite  values .  We  attain  this  by  assuming  in  the  fluid  so-called 
Helmholtz  discontinuity  surfaces  which  extend  from  the  plate  edges.   By 
this  expression  one  understands  vortex  sheets  consisting  of  free  vortices 
and  obeying  Helmholtz'  laws.   In  contrast  to  this  is  the  plate,  which 
also  is  a  discontinuity  surface  but  no  longer  obeys  Helmholtz'  laws.   It 
forms  a  rigid  surface  and  exerts,  therefore,  pressures  on  the  flow. 

The  plane  in  which  the  flow  takes  place  will  henceforth  be  designated 
as  complex  z-plane.   The  system  of  coordinates  is  selected  so  that  the 
origin  of  the  coordinates  coincides  with  the  center  of  the  plate  and  the 
plate  lies  along  an  imaginary  axis.   The  width  of  the  plate  is  made  equal 
to  2b .   The  movement  of  the  plate  is  to  take  place  from  the  left  to  the 
right,  or,  which  means  the  same,  the  flow  is  to  approach  the  plate  from 
the  positive  side.   The  constant  velocity  at  infinity  -  which,  with  oppo- 
site sign,  may  also  be  interpreted  as  speed  of  the  plate  traveling  in  a 
fluid  at  rest  -  is  to  be  -v^      (fig-  1^  left).  We  conformally  map  the 
z-plane  onto  a  ^-plane,  by  means  of  the  transformation 


/ 


2  +  b2  (1) 


with  the  plate,  which  represents,  of  course,  a  piece  of  the  imaginary 
axis  in  the  z-plane,  becoming  a  piece  of  the  real  axis  in  the  new  plane 
(fig.  1,  right).   The  edges  z   =   +ib  shift  into  the  zero  point  ^  =  0; 
the  point  at  z  =  0  transforms  to  the  points  ±b  on  the  real  ^-axis. 
Furthermore,  in  this  transformation  the  region  at  infinity  in  the  z-plane 
goes  over,  without  change,  into  that  at  infinity  in  the  ^-plane.   There- 
fore the  velocities  at  infinity  in  both  planes  are  equal,  thus 
v(z  =  Qo)  =  v(^  =  00)  =  -v^.   As  images  of  the  discontinuity  surfaces  we 
obtain  in  the  new  plane  again  images  which  start  from  the  point   ^  =  0. 
The  new  vortex  sheets  lie  symmetrically  to  the  real  axis  whereby  the 
latter  becomes  a  streamline  so  that  the  image  of  the  plate  also  remains 
a  streamline-'-.   The  flow  in  the  ^-plane  corresponding  to  the  plate  flow 
is  composed  of  the  parallel  flow  and  a  flow  caused  by  the  free  vortices 
of  the  discontinuity  surfaces.  We  are  now  going  to  calculate  its 
velocity  field. 

For  the  case  where  the  vortex  sheets  in  the  z-plane  no  longer 
lie  symmetrical  to  the  real  axis,  one  must  choose  instead  of  (l)  the 
transformation  into  a  circle  since  the  real  axis  in  the  ^-plane  then  is 
no  longer  a  streamline.  For  our  case  the  selected  transformation  is 
preferable  to  that  into  the  circle  because  the  graphical  calculations 
become  considerably  simpler. 
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We  denote  by   ^-p  a  variable  point  of  the  vortex  sheet  in  the  upper 

half  plane  and  coordinate  to  it  a  circulation  dP.  The  corresponding  vortex 

point  of  the  mirrored  surface  is  then  given  by   ^-p,   and  has  the  same 

vortex  strength,  but  with  opposite  sign,  thus  -df.  Both  vortex  points 
produce  at  an  arbitrary  point  ^  of  the  plane  a  velocity  the  conjugate 
complex  value  of  vhich  is  given  by 


dv^(0  = 


dT 


2.1^^ -^F     ^_g 


If  one  puts  dr=  7(^p)d^p  wherein  /(^p)  represents  the  distribu- 
tion density  at  a  point   ^p,  and  d^„  a  line  element  of  the  surface,  one 

obtains  by  integration  over  all  points  of  the  surfaces  the  velocity 
induced  by  them 


(^p'   signifies  the  temiinal  point  of  the  vortex  sheet.)  The  parallel 

flow  has  the  velocity  -v^.  For  the  total  velocity  of  the  superimposed 

flow  one  may  now  write 


=F        (2) 
Between  the  velocities  of  both  planes  there  exists  the  relationship 

H^  =  HU^  (3) 

dz 


from  which  one  obtains  the  velocities  for  the  z-plane.   The  condition 
that  the  velocity  should  assume  finite  values  at  the  plate  edges  remains 
to  be  satisfied.   As  one  can  see  very  quickly,  this  requirement  is  ful- 
filled when  the  expression  (2)  for   C  =  0  disappears,  because  the  magni- 
fication ratio  d^/dz  becomes  infinite  for  •  ^  =  0  or  z  =  tib .   One 
obtains  therewith  the  condition^ 

The  part  stemining  from  V]_(^)  becomes  real  for   ^  =  0. 
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which  must  always  be  satisfied.   It  is^  however,  not  sufficient  for  deter- 
mination of  distribution  density  and  shape  of  the  discontinuity  surface. 

As  long  as  the  vortices  are  still  very  smaH  compared  to  the  plate, 
the  plate  may  be  regarded  as  infinitely  wide  in  comparison  with  the  vor- 
tices.  For  this  case,  the  lack  of  a  fixed  comparative  length  is  the 
reason  that  the  form  of  the  discontinuity  surface  is  independent  of  the 
magnitude  of  the  vortex  configuration.   The  development  of  the  discon- 
tinuity surface  consists,  therefore,  initially  in  a  similar  magnification. 
This  similar  magnification  can  be  fulfilled  only  in  the  case  of  a  certain 
shape  of  the  discontinuity  surface,  and  of  a  certain  distribution  of  the 
vortices  over  it.   On  the  basis  of  this  condition,  we  can  therefore  cal- 
culate the  form  and  vortex  distribution  for  the  beginning  of  the  motion. 
Starting  from  this  initial  condition  we  can  then  find,  by  further  obser- 
vation of  the  variation  with  time  of  the  change  in  form,  'also  the  shapes 
and  vortex  distributions  for  the  later  times  when  the  vortices  are  no 
longer  small  compared  to  the  width  of  the  plate . 

Before  performing  the  limiting  process  to  the  plate  of  infinite 
width,  we  choose  a  new  system  of  coordinates   z'   which  arises  from  the 
former  coordinate  system  by  parallel  displacement,  and  the  zero  point  of 
which  is  shifted  to  the  upper  edge  (z  =  +ib).  We  write 


z  =  z'  +  ib  (5) 

and  replace,  in  addition,   z   by 

z'  =  5z  or  Z  =  ^  z'  (6) 

b  H 


where  H  represents  a  scale  unit.   Accordingly,  equation  (l)  assumes 
the  form 


^  =  l/fSz')^  +  2iHZ  (7) 


The  limiting  process  from  the  plate  of  finite  width  to  that  of  infinite 
width  (b  — >  00)  yields  a  new  traas  format  ion 

I t^ 

^  =  \J2±EL     or  Z  =  -  i  |jj  (8) 
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This  is  the  transformation  of  the  plane  with  a  slit  along  the  negative 
imaginary  axis  onto  the  upper  half  plane  of  the  ^-plane .   If  in  the 
^-plane  v^^  continues  denoting  the  free-stream  velocity,  the  interpre- 
tation of  v^  as  traveling  or,  respectively,  free-stream  velocity,  is 

lost  in  the  Z-plane  since  in  this  plane  the  velocity  at  infinity 
tends  ^0.   The  quantity  v^^  appears  in  this  plane  as  the  velocity 

which  prevails  in  the  case  of  vortex- free  flow  at  the  point  Z  = iH. 

2 

All  other  considerations  carried  out  so  far  concerning  the  plate  of 
finite  width  may  he  directly  transferred  to  the  plate  of  infinite  width. 


5.  THE  INFINITELY  WIDE  PLATE 
(a)  Laws  of  Similitude 


By  the  limiting  process  from  a  finitely  wide  to  an  infinitely  wide 
plate  which  was  just  completed,  a  flow  type  was  obtained  for  which  the 
successive  flow  patterns  are  similar  and,  accordingly,  the  circulation 
of  the  vortex  and  the  distribution  density  on  the  vortex  sheet  are 
similarly  enlarged  for  this  growth.   We  shall  briefly  derive  here  the 
laws  which  characterize  this  behavior,   (As  was  remarked  in  the 
Introduction,  the  laws  of  similitude  could  be  immediately  derived  from 
the  quoted  treatise  of  Prandtl.) 

We  observe  the  distance   Z  of  two  points  remaining  in  similar 
position  during  the  similar  enlargement  of  the  discontinuity  surface  and 
assume  it  to  be  at  the  time  t-,   of  the  amount   Z-|_  and  at  the  time  to 

of  the  amount   Z2.   We  state  as  the  time  law  for  the  enlargement 

Z2   \t2J 

where     A     is   a  number  yet  to  be  determined.      The   law  for  the   circulation 
r     about   similarly  situated  regions  may  be   obtained  as   follows.      Since 
the   free-stream  velocity  in  the    ^-plane   is   invariably     Vqo,    the   same  veloc- 
ity must  prevail  also   for  the  similar  enlargement   at   similarly   located 
pointG5.      This   velocity  is   composed  of  the   free-stream  velocity     v^o     and 

of   the   fields   of  the   individual  vortex  elements.      The   influence   of   a 

In  the  Z-plane   this   is  not  the   case,    since  here    (for  vortex-free 


ii: 
2 

during  the  similar  enlargement,  thus  does  not  remain  in  similar  position. 


flow)  the  velocity  is  constant  at  the  point  'Z  =-—   iH  which  is  fixed 
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vortex  element  AF^  at  the  distance   Zf   from  the  point  considered^  on 
the  velocity  is  Ar/2nlj-.      If  during  the  gro-wth  AF  and   Zf   are  modified, 
AF/Zf,  due  to  the  constant  velocity,  also  must  be  constant.  For  the 
circulation   F  about  similarly  situated  regions  of  the  ^-plane,  therefore, 


F2  z^ 


2 


is  valid.   In  the  Z-plane  the  corresponding  distances  axe 


^1  _  hi"^' 


2/ 


Z2         yZ^ 
Thus     —    becomes  in  the   Z-plane 

r2 


^,t  © 


If  we  designate  by  V  the  velocities  which  correspond  to  the  law  of 
similitude,  we  obtain  for  the  velocity  V(z)  in  the  Z-plane  which  is  pro- 
portional to   r/z 


Vi_^iAi_    1^     /tiN-^/2 


Fg/Zg    \j  Z-]_   \t^i 


(11) 


If  we  replace  in  the  ratios  above  the  time  t-|_  by  the  time   1  and  the 
time  tp  by  the  time  t,  the  laws  of  similitude  read 

z  =  z^t^   F=  r^t^s  V  =  Vit"^/2  (12) 
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wherein  the  number  A  remains  to  be  determined.  For  calculation  of  the 
number  A  we  form  the  velocity  of  growth  of  the  distance  I 


dZ    ,  ,  .  A-1 

=   6  T  At 

dt    -^ 


Since  this  velocity,  like  all  velocities,  must  vary  proportionally 

-A/2 
to  t~   '    ,    there  results 

A-l  =  -^ 


thus 


A^2 


4 


Therewith  the  laws  of  similitude  read 

z  =  Zit2/5      r  =  rit^/5      V  =  Vit-1/5  (15) 

(b)  Solution  for  the  Vortex  Core 

Our  task  now  consists  in  finding  a  shape  for  the  spiral  which  satisfies 
these  laws.   They  are  the  same  laws  Kaden  (ref.  5)  obtained  as  a  result. 
Thus  one  deals  here,  too,  with  a  similar  problem.   In  particular,  both  prob- 
lems become  identical  for  the  inner  part  of  the  vortex  where  in  both  cases 
a  discontinuity  surface  consisting  of  free  vortices  is  to  be  rolled  up  into 
a  spiral  vortex.   It  is  therefore  sufficient  if  we  refer  to  Kaden 's  paper 
and  here  only  briefly  mention  the  results. 

Let  R  and  cp  be  polar  coordinates  of  a  system  of  coordinates  the 
zero  point  of  which  coincides  with  the  center  of  the  spiral.   For  form  and 
circulation   F  of  the  vortex  core  bounded  by  a  circle  of  the  radius  R, 
Kaden  obtained 


^=m^/3,     r  =  2xv/S  (11+) 
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The  circumferential  velocity^  for  radius  R  is 


^^  =  ^  =  ^  (15) 


Therein  t  is  the  time,  counted  from  the  beginning  of  the  motion.   The 
quantity  X  is  a  constant,  still  unknown  for  the  time  being.   It  has  the 
dimension  velocity  times  square  root  of  a  length.   Thus  we  may  also  intro- 
duce for  it  the  dimensionless  constant 


k  =  -^  (16) 


(c)  Solution  for  the  Outer  Loops  of  the  Spiral 

As  was  shown  before,  the  behavior  of  the  spiral  is  characterized  in 

/xt\2/3 
the  interior,  that  is,  for  large  angles  9,  by  R  =  —    .  For  small 

angles  cp,  in  contrast,  especially  where  the  discontinuity  surface 
adjoins  the  edge,  considerable  deviations  from  this  form  occur.   This 
transitional  region  can  be  found  as  follows.   We  visualize  the  plate  as 
lying  parallel  to  the  straight  line  cp  =  0;  this  assumption  is  insignifi- 
cant for  practical  purposes  since  the  spiral  windings  approach,  in  the 
inward  direction,  a  circular  form.   We  isolate,  furthermore,  at  one  point 
the  inner  part  of  the  spiral,  in  which  the  form  is  prescribed  with  suffi- 

/xt^2/5 
cient  accuracy  by  R  =  —  I   .   For  calculation  of  the  velocity  field 

Vjtqpy 

outside  of  the  spiral  core,  we  visualize  the  latter  as  replaced  by  an 

isolated  vortex  at  the  center  of  gravity  of  the  circulation.   This  is 

directly  permissible  since  almost  circular  symmetry  prevails  in  the 

interior  of  the  spiral.   The  dissymmetry  caused  by  the  separation  point 

brings  it  about  that  the  center  of  gravity  of  the  circulation  does  not 

coincide  with  the  center  of  the  spiral.  For  the  calculations  indicated 

later  on,  the  separation  point  is  placed  at  the  point  cp  =  2.5fl:«   The 

difference  between  the  coordinates  of  the  spiral  center  (a,  h)  and  of  the 

center  of  gravity  (a;i_,  h^)   may  be  calculated,  similarly  as  by  Kaden, 

(ref.  3)  to  be  approximately 
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a^  -  a  =  0.006  R^   and  h^  -  h  =  O.O5I  R^  (17) 

if  R.   is  the  radius  of  the  core  of  the  spiral. 

The  outer  loop  which  we  now  want  to  calculate,  extends  from  the 
separation  point  (cp  =  2.5rt)  to  the  plate  edge  (cp  =  cp*).   Its  shape  and 
vortex  distribution  7  as  well  as  the  position  and  circulation  F^     of 

the  vortex  core  we  assume,  at  first,  arbitrarily,  with  a  factor  of  simili- 
tude  X  common  to  both  quantities  and  at  present  not  yet  determined, 
(equations  (ik)   and  (I5)),  and  make  therewith  the  transition  into  the 
^-plane.   Therein  7  is  transformed  into 


7'  =  7  ^ 


whereas  the  circulation  rj_     remains  unchanged.   In  this  mapping  plane 

(^-plane)  the  velocities  stemming  from  the  vortex  may  be  calculated  as 
the  field  of  the  vortex  core  and  the  distribution  on  the  outer  winding, 
likewise  their  mirrored  images.   By  superimposltion  of  the  undisturbed 
velocity  v^  there  results  the  velocity  field  v'   in  the  ^-plane. 

First,  one  ascertains  the  velocity  at  the  point   ^  =  0.  From  the  condi- 
tion that  this  velocity  must  be  zero  (equation  (k))   results  factor  X, 
still  undetermined  at  first,  for  the  circulation  f^   and  the  distribu- 
tion density  7. 

In  the  conformal  mapping  onto  the  Z-plane  the  velocities  v'   of 
the  ^-plane  are  transformed  into  the  velocities  v  of  the  Z-plane, 
according  to  the  relationship 

-   -,  d^ 

V  =  V  — 

dZ 


where  v  and  v   signify  the  conjugate  values  of  the  velocities.   The 
velocity  of  the  individual  vortex  concentrated  at  the  core  is  obtained 
by  omission  of  the  field  of  this  vortex.   It  must,  however,  be  noted 
that  in  the  transition  from  the  ^-plane  to  the  Z-plane  the  field  of  this 
vortex  is  deformed;  hence  an  additional  term  appears  in  the  conversion 
of  the  velocities.   In  the  ^-plane  the  potential  of  the  core  vortex  to  be 
omitted  is 


1 
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*k(?)  —!-(?-  to) 

In  the  Z-plane  the  potential  of  the   vortex  to  be   omitted  is 

23ti         V  ^0  / 


2fli  V''         ^0  2:ti         V  0 


In  order  to  calculate  the  velocity  for  the  Z-plane  we  must^  therefore, 
subtract,  besides  the  potential  ^y-i^)  ^   also  the  additional  potential 

r. 

— —   In  K  +  t^   \      to  which  corresponds  a  velocity 
2jti    V     'J/ 


v' 

additional 


in  the  ^ -plane  or,  respectively. 


^i 


2iti  U  +    C 


"1  rlr' 


2n±  K  +   U    dZ 


additional 


in  the  Z-plane . 

Once  one  has  calculated  the  velocities  for  individual  points  of  the 
discontinuity  surface  one  resolves  them,  most  advantageously,  into  normal 
and  tangential  velocities  v^^  and  v^,  and  plots  the  latter  as  functions 

of  the  arc  length  s .   By  combination  of  the  function  values  one  obtains 
the  velocities  pertaining  to  each  point  of  the  vortex  sheet;  from  them 
results  the  motion  of  the  assiomed  vortex  sheet. 

A  fluid  particle  at  the  point  P  of  the  spiral  having  the  velocity 
V  (fig.  2)  moves  by  the  distance  PP-i  =  v  in  the  unit  time.   Due  to  the 


12 
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similar  enlargement,  the  point  P  is  tranformed  during  this  time  into 
the  similarly  located  point  P2« 

The  velocity  at  which  the  similarly  situated  points  are  displaced 
is  assumed  to  be  V  so  that  the  distance  is  PP2  =  V.   To  have  P-^ 

come  to  lie  on  the  similarly  enlarged  spiral,  the  normal  components 
V   and  Vj^  of  the  two  velocities  v  and  V  must  be  equal  or 


AVn  =  v^  -  Vn  =  0  (19) 


We  observe  the  point  P  on  its  path  to  the  similarly  situated 
point  Pp  and  find  that  this  path  does  not  represent  the  motion  of  a 
fluid  particle.  Rather,  the  fluid  flows,  with  the  velocity  v  -  V, 
through  the  point  moving  from  P  to  P2. 

If  V1J1  and  Vrp  signify  the  tangential  components  of  the  velocities 

V  and  V,  the  circulation  in  flowing  with  the  fluid  into  the  core  inside 
of  P  or  P2^  respectively,  during  the  time  At  is 

Af  =  (vt  -  V^)7At 


Hence  there  results  for  the  region  inside  of  P,  which  is  being  similarly 
enlarged 

—  =  (vrp  -  V^)/ 
ot 


According  to  the  laws  of  similitude  (13)^  one  must  have  for  a  region 
being  similarly  enlarged 

^r  ^  r  1  t-2/5  =  _L 
St    ^  5       5t 


The  two  values  of  dP/dt  must  agree  in  order  to  have  the  similitude 
satisfied,  therefore  the  equation 
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A  =  (Vrn 

dt      ^ 


Vt)7 


3t 


=  0 


(20) 


must  apply.   Since  the  distribution  density  is 


\ 


we  may  write  instead  of  (20)  also 


ht  OS 


(Vrp  -  Vrp)7 


^  =  0 

3t 


(21) 


With  the  aid  of  the  equations  (19)  and  (2l)  it  is  now  possible  to  find 
the  shape  and  vortex  distribution  of  the  outer  loop. 

We  shall  perform  the  calculation  of  these  quantities  for  a  time  t 


t  =  «  h5/2 
X 


(20)' 


for  which  the  equation  (ik)   of  the  spiral  core  assumes  the  simplified 
form 


R 


H 


,2/5 


(21)' 


9' 


We  select  as  the  initial  form  for  the  entire  spiral  the  form  R  =  H/cp^/^ 
and  place  it  as  a  first  try  so  that  It  joins  the  plate  at  the  point 
cp  =  cp-x-  =  jt/2  (fig.  5).   The  break  originated  at  the  edge  is  merely  a  flow 
which^  as  one  readily  understands,  disappears  at  the  next  moment.   Except 
for  the  factor  x^  the  circulation 


2XJR 


^ 


(22) 




The  duplication  of  the  equation  numbers  20  and  21  follows  that  of 
the  original  Germain  docujnent. 
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is  known  to  us.   For  the  core  which  we  shall  assume  to  begin  at  cp  =  2.5it, 
we  have 

R  =  R   =  1_^ 


(2.5rt) 


Thus  r.   becomes 


r,  =  2x,.  / 5 (25) 


V(2.5«)2/5 


The  distribution  density  is 


7 


-  Sr  _  ar  dR  -  V  R       .  ,  X  ^21,) 


Bs        Sr  Ss 


\/9Ah5  +  r5      \/icp2/3\|9/i^  +  cp- 


In  figure  4  this  distribution  density  7  is  plotted  as  a  function 
of  the  arc  length  s.   The  conforraal  mapping  (8),  in  combination  with  the 
reflection  at  the  |-axis,  leads  to  the  double  spiral  indicated  in  figure  5> 
The  distribution  density  in  this  plane  (image  plane)  is  represented  in 
figure  6.  When  we  calculate  the  velocity  at  the  zero  pointy  we  obtain, 
on  the  basis  of  equation  {h^ y    8-  condition  for  the  still  undetermined 
value  X.   There  results 


<-L  =  2.36\/Hv„  or  k-[_  =  2.56 


As  an  example  for  the  determination  of  the  velocity  field  we  shall 
here  briefly  reproduce  the  calculation  of  the  velocity  at  the  center  of 
gravity.   In  the  ^-plane  an  element  of  the  discontinuity  surface  ds 
induces  at  the  center  of  gravity  the  velocity  components 

.        /(s '  )  ,  ,        /(s '  )  , 

dvt        =  sin  \|/  ds    ,      dv^,        =  cos  ^  ds 

^1  2jtRs'  '1  2jtRs' 
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when  Rg '   denotes  the  distance  of  the  element  from  the  center  of  gravity 

and  ijr  the  angle  of  the  radius  vector  from  the  center  of  gravity  to  the 
element  with  the  ^-axis  (fig.  5)«   By  integration  over  s'   from  the  zero 
point  to  the  core  boundary  (cp  =  2.5jt,  s'  =  2.62  H)  and  over  the  corre- 
sponding curve  (s')  of  the  image  one  obtains  vt  '   and  "v^  '  .   In  fig- 
ure 7  the  values  dv^  '/ds'   and  dv^  '/ds'   are  plotted  against  s'   or 
s '■  respectively.   The  quantities  ve  '   and  v^  '   are  obtained  by  cir- 
cumscribing; with  the  aid  of  a  planimeter^  the  cross-hatched  areas 
bounded  by  the  curves  and  the  abscissa  axis 


These  values  represent  the  influence  of  the  outer  loop  of  the  spiral. 
The  image  of  the  vortex  core  causes  the  velocities 


Ti 


=  0.22v_     v„   =  0 


^2    2«2h'     ^   ~      ^2 


if  h'   signifies  the  distance  of  the  center  of  gravity  of  the  core  from 
the  |-axis.   The  vortex  core  itself  does  not  contribute  to  the  velocity 
of  its  own  center  of  gravity  in  the  ^ -plane.   However^  it  must  be  noted 
that,  according  to  the  explanations  on  pp.  10  and  11,  in  the  conversion 
to  the  Z-plane  a  term  stemming  from  the  vortex  core 


vt^'  =  0.11  v_     Vr,,  '  =  0.11  v„ 


^13  -  u.xi  v„     v^^ 


appears , 

For  the  resultant  velocity  one   obtains 


^1  '    =  -  ^co  +  "^li  +  v^2  ^  "^^3   "  -O.i^Ov^ 


Vti     =  v^^  +  v^^  +  v^^  =  0.11v„ 
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whence  there  results,  by  conversion  to  the  Z-plane, 

Vx  =  H  —^ f-   =  -  0.29v<„ 

r   +  ^ 


^y  =  H-^- -^=  -  0.17v„ 

r   +  r\ 


According  to  'the  Law  of  similitude, 


Vx  =  -  0-57v<„     Vy  =  0 


would  have  to  be  valid.   The  center  of  gravity  thus  shows  for  the 
assumptions  made,  relatively  to  its  required  motion,  a  wrong  motion 
obliquely  downward,  toward  the  plate,  with  the  velocity  components 


^^x  =  ^x  =  V^  =  0-08v^ 


^V  =  ^y  -  Vy  =  -  0.l7v„ 


In  the  same  manner  the  velocities  for  the  points  of  the  outer  winding  are 
calculated.  Of  course,  one  has  to  add  a  term  which  represents  the  influ- 
ence of  the  core;  on  the  other  hand,  the  additional  term  V:^  does  not 

appear  here.   We  find  the  calculated  normal  and  tangential  components 
plotted  in  figures  8a  and  8c.   As  can  be  seen  from  figures  8b  and  8d, 
equations  (I9)  and  (21)  are  not  satisfied.   In  other  words,  the  shape  and 
vortex  distribution  assumed  deviate  from  the  actual  solution.   We  now 
modify  the  position  of  the  core  e.g.  in  the  sense  of  Av^^  and  Av„  and 

the  shape  and  vortex  distribution  of  the  outer  loop  in  the  sense  of  the 

differences  Av„   and  A  -^.      The  fonn  of  the  spiral  core  remains 
"^        dt 


'Xt\  2/5  ^   H 


unchangedly  R  =  (-i_i     =  .   We  now  repeat  the  calculation  procedure 

H  cp5/2 
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carried  out  so  far.   After  several  steps,  the  first  two  of  which  are 
represented  in  figures  9  "to  12,  one  arrives  finally  at  a  solution  that 
satisfies  the  laws  of  similitude;  it  is  plotted  in  figure  13 • 

In  figure  li<-  is  the  radius  of  the  spiral;  in  figure  I5  its  distri- 
bution density  is  plotted  as  functions  of  the  center  angle  and  of  the 
arc  length  s,  respectively. 

For  more  convenient  further  use,  the  most  important  constants 
necessary  for  characterization  of  the  spiral  vortex  have  heen  briefly 
compiled  below.   For  the  solution  indicated  in  figures  15  to  I5,  the 
calculation  yields  the  time 

t  =  1.42  i 


the  constant 

X  =  2.22H^/^      k  =  2.22 
00 

the  coordinates  of  the  spiral  center 

a  =  -0.53H     h  =  -0.22H 
the  total  circulation  of  the  spiral  vortex 

r  =  lj-.20Hv„ 


For  the   assiomed  core  boundary  at     cp  =  2.5rt      the   radius   of  the   core 
is 


R^   =  O.267H 


the   circulation  of  the   core 


Ti   =   2.3Hv„ 
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For   an  arbitrary  time     t     one   obtains 

a  =   -0.i+Ht2/3  h  =   -0.l65Ht2/5 


R  = 


r  =  3.7Hv<„t 


=    ^A.    TT+2/5 


1/3 


a 


cp2/5  cp2/3 


Ri   =  0.2Ht^/5  =  0.5a 


r.    =  2Hv^t^/5  ^  5.i5v/I 


i4-.  THE  PLATE  OF  FINITE  WIDTH 


With  the  solution  found  for  the  plate  of  infinite  width,  we  know 
the  initial  flow  conditions  about  a  plate  of  finite  width,  where  the 
vortices  are  still  small  compared  to  the  plate  width.   We  know  that  for 
this  initial  state  the  successive  vortex  images  originate  from  one 
another  by  similar  enlargement.   However,  as  soon  as  the  vortices  assume 
a  magnitude  which  is  no  longer  negligible  compared  to  the  plate  width, 
the  presupposition  for  a  similar  growth  no  longer  holds  true.   Yet, 
starting  from  the  form  and  distribution  density  we  fovmd  for  the  still 
very  small  vortices,  we  can  calculate  the  variation  of  this  form  and 

distribution  density,  using  the  values  of  Av^^  [(19)^  and  A  —  (2l) 

which  now,  with  growing  vortices,  more  and  more  deviate  from  0.   We  only 
must  make  use  of  the  transformation  by  means  of  (l)  instead  of  the  trans- 
formation by  means  of  (8).   If  we  start  from  a  known  form  and  distribu- 
tion density  at  a  time  t-i  ,  we  obtain  the  deviation  of  the  form  and  dis- 
tribution density  from  those  which  would  result  for  similar  enlargement 
according  to  the  time  laws  (13) ,    at  the  time  to  as  being 


t2  n  tg 


by 


An  =      /  AVj^dt      and     Aj  =  A  — ^  dt 


^^1  -^^1 


bt 


I 
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The  graphical  determination  of  the  required  quantities  is  greatly 
hampered  by  the  fact  that  the  vortices  are  at  the  beginning  very  small 
and  later  very  large.   We  can  avoid  this  difficulty  by  visualizing  the 
figure^  at  every  instant,  enlarged  or  reduced  in  such  a  manner  that  the 
vortex  would  always  remain  the  same  if  it  would  keep  on  growing  according 
to  the  laws  of  similitude  that  are  valid  for  the  beginning.   The  figure 
in  actual  size  we  denote  as  figure  I  and  the  quantities  valid  for  it  by 
the  subscript  I,  the  enlarged  figure  as  figure  II  and  the  corresponding 
quantities  by  the  subscript  II.   At  the  initial  stage,  the  enlarged  vor- 
tex is  to  agree  precisely  with  the  one  calculated  for  the  plate  of 
infinite  width  in  the  preceding  section.   Let  us  call  the  ratio  of  mag- 
nification in  every  case   e.   If  we  reduce  in  the  magnified  figure  all 
velocities  simultaneously  at  suitable  points  in  the  proportion  l/v/e", 

the  circulations   F  are  magnified  in  the  proportion  \fe.      The  distribu- 
tion densities  7  are  reduced  in  the  proportion  l/u^- 

For  the  potential  flow  about  the  plate  of  the  width  2b  in  the 
plane  I  there  results  in  the  neighborhood  of  the  edge  the  velocity 


TI  =  -pj=l=  (25) 


^^ 


For  the  plate  of  Infinite  length  treated  in  section  3  "the  corresponding 
velocity  was 


(26) 


wherein  y'  =  b  -  y  signifies  the  distance  from  the  edge  of  the  plate. 
In  order  to  obtain  for  the  beginning  of  the  motion  the  same  conditions 
as  for  the  flow  about  the  plate  of  infinite  width  treated  before,  we 
select  the  width  of  the  plate   2b  in  such  a  manner  that 


b  =  H 

For  the  above  stipiilation,  exactly  the  velocity  V;„  prevails  in  the 

magnified  figure,  while   e  is  still  very  large,  at  a  point  at  a  distance 
of  b/2  from  the  edge.   In  the  case  of  the  plate  of  infinite  length 
treated  before,  this  point  lay  at  a  distance  H/2  from  the  edge.   Since 
we  equated  b  =  H,  we  obtain  in  the  magnified  figure  precisely  the  flow 
treated  in  section  3. 
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According  to  equation  (1^4-)^  for  the  growing  vortex  core  of  the 
plate  of  infinite  width 


Since  at  the  initial  stage  the  flow  about  the  edge  of  the  plate  of  the 
finite  width  2b  =  2H  coincides  with  the  flow  about  the  edge  of  the 
plate  of  infinite  width^  this  formula  applies  also  to  the  initial  stage 
in  the  case  of  a  plate  of  finite  width.   By  means  of  the  magnification 
by  the  factor  e,  this  initial  vortex  is  to  be  transformed  into  the  flow 

TT       >, 

treated  for  which     R   =  — ^j-  =  j—.      Thus  necessarily 


q)2/3       cp2/3' 


Re  =   e 


Xt\ 


2/5 


«^/  cp2/3 


Hence   there   results 


.  =  (1)2/5  M-2/5  (^7) 


Since  in  the  plane  I^  at  the  beginning  of  the  motion,  the  circulations 

1/5 
Fj     grow,  according  to  (l3),   with  t  '  ,  we  obtain  in  the  plane  II,  where 

they  appear  magnified  in  the  proportion  ye,  constant  circulations  and 
constant  distribution  densities. 

At  the  beginning  of  the  motion  e  =  <».   The  plate  edges  lie  in  the 
plane  II  at  an  infinite  distance  from  one  another  so  that  we  actually 
have  the  case  of  the  plate  of  infinite  width.   In  time,  however,   e 
attains  smaller  finite  values,  and  we  obtain  in  the  plane  II  also  a  plate 
of  finite  width.   Therewith  the  velocities  become  different,  and  we 
obtain  deviations  from  the  similar  magnification.   In  order  to  calculate 
the  velocities  in  the  plane  II,  we  transform  this   Zjj-plane  onto  a 

^-plane  which  we  shall  characterize  by  the  subscript  III,  by  means  of  the 
function 


te  =  lP^- (28) 
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For  very  large   e  this  mapping  is  transformed  into  the  infinite  plate 
according  to  (8)  in  the  neighborhood  of  the  edge.  For  finite  e, 
however,  other  shapes  result  for  the  mapped  vortices,  and  that  is  the 
very  reason  which  causes  the  modifications  of  the  velocities  and  there- 
with the  deviations  from  similar  vortex  growth. 


As  soon  as,  due  to  this  deviation,  the  form  and  distribution  density 
of  the  vortices  have  changed  with  respect  to  the  similar  growth,  this 
also  contributes  to  the  variation  of  the  velocities.   Since,  however, 
form  and  distribution  density  change,  at  first,  only  very  slowly  compared 
to  the  similar  growth,  one  may  in  the  plane  II  assume  the  form  and  distri- 
bution density  as  constant  in  turn  through  a  large  time  interval,  and 
need  consider  in  this  time  interval  only  the  modification  of  the  trans- 
formation in  the  plane  III,  due  to  the  modification  of  the  value   e. 

For  such  a  time  interval  t-^  to  t2  (the  first  starts  with  t  =  0), 

one  calculates  for  several  intermediate  times  the  normal  and  tangential 
components  Vn^-  and  v-^    of  the  velocities  of  the  plane  II  with  the 

aid  of  the  transformation  onto  the  plane  III  in  the  same  manner  as  in  the 
case  of  the  plate  of  infinite  length.   One  forms  furthermore  the  differ- 
ences with  respect  to  the  velocities  of  the  similar  magnification  Vr 


and  V+ 


^11 


'nil 
and  obtains  then  by  transfer  to  the  plane  I  the  values  (l9) 


Av 


^I 


t-'-ii '  ^f-ii 


=   €Av 


-  V, 


n 


11/ 


(29) 


and  corresponding  to  (21) 

3t    Ss-j- 

(-Ti  -  Vtj)7i 

71 

3t 

=  -e  

Ssj 

I 

[    ^11      ^13 

)/^-/ 

3t 


5/2 

If  we,   finally,   replace   in  the   last  term     t     by   (b/e)    '      Jt/x  (equation 
(27)),  we  obtain 


B7  . 

A— 1  =  -e' 
Bt 


=11 


■II 


■^Ilj  ^11  ^ 


X7ii 


3jtb 


3/2 


(30) 
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By  graphical  integration  then  results  the  displacement  An  of  the  points 
of  the  vortex  sheet  at  right  angles  to  it,  and  the  modification  of  the 
distribution  density  compared  to  the  similar  magnification  in  the 
plane  I  as 


Z\ni 

ti 

^  Av^^dt 

A7,   = 

Jti 

A  — i  dt 

(31) 


(52) 


and  the  modification  of  form  and  distribution  density  in  the  plane  II 
as 


Anjii  =  egAnj  (53) 


^''11  =  ^  %  (54) 


Therein  £2  is  the  value  of  the  magnification  ratio   e  at  the  time  t2. 

Due  to  the  finite  magnitude  of  the  vortices  compared  with  the  plate 

^7ti 
width,  Av„   ^   0  and  A  f   0:  also,  the  condition  (k) , 

''ll  at 

that  in  the  ^ -plane  at  the  zero  point  the  velocity  must  be  zero,  will 
no  longer  be  satisfied.   In  the  ascertainment  of  the  vortex  for  the  plate 
of  infinite  width  we  have  been  able  to  fulfill  this  condition  by  suitable 
definition  of  an  as  yet  undetermined  factor  for  the  circulation.   Due 
to  this  condition  we  found  the  quantity   X  or   k,  respectively.  For 
the  further  development  of  the  vortex,  fonii  and  vortex  strength  and 
their  variation  with  time  are  fixed.   Only  the  strength  of  the  vortices 
shed  at  the  plate  edge  is  still  undetermined  since  we  can  here  not  form 
the  differential  quotient   S/Ss  occurring  in  (50).   We  must  select  it 
in  such  a  manner  that  equation  (h)    is  satisfied.   We  obtain  therefore  an 
additional  modification  of  the  distribution  density  starting  from  the 
plate  edge  which  gradually  is  carried  into  the  vortex  by  the  flow.   It 
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is  true  that  it  was  shown  in  the  quantitative  calculation  that  the  devia- 
tions from  the  condition  (k)   are  extremely  small,  because  the  variation 
of  the  conformal  representation  [(l)]  with  e  results  in  a  positive 
velocity,  whereas  the  modification  of  the  form  and  distribution  density 
according  to  (33)  and.  (3^)  results  in  a  negative  velocity  at  the 
point   ^  =  0  and  the  two  almost  cancel  one  another. 

The  calculation  was  carried  out  for  the  intervals   e  =  00  to  e  =  3> 
e  =  3  to   6=2,  aind  e  =  2  to  e  =  1.   The  results  are  compiled  in  the 
figures  16  to  I9.   True  to  expectation,  the  circulation  F   increases 
more  slowly  with  time  than  it  does  according  to  the  solution  for  the 
plate  of  infinite  width  (initial  state,  fig.  19) •   In  the  final  state  it 
would  perhaps  approach  a  constant  value  which  corresponds  to  a  steady 
state  of  flow.   However,  according  to  experience  the  symmetrical  vortex 
configuration  becomes  unstable  from  a  certain  magnitude  onward,  so  that 
this  steady  state  is  not  attained. 


5.  SIMVIARY 


The  flow  about  the  plate  of  infinite  width  may  be  represented  as  a 
potential  flow  with  discontinuity  surfaces  which  extend  from  the  plate 
edges.   For  prescribed  form  and  vortex  distribution  of  the  discontinuity 
surfaces,  the  velocity  field  may  be  calculated  by  means  of  a  conformal 
representation.   One  condition  is  that  the  velocity  at  the  plate  edges 
must  be  finite.   However,  it  is  not  sirfficient  for  determination  of  the 
form  and  vortex  distribution  of  the  surface.   However,  on  the  basis  of 
a  similitude  requirement  one  succeeds  in  finding  a  solution  of  this 
problem  for  the  plate  of  infinite  width  which  is  correct  for  the  very 
beginning  of  the  motion  of  the  fluid.   Starting  from  this  solution,  the 
further  development  of  the  vortex  distribution  and  shape  of  the  surface 
are  observed  in  the  case  of  a  plate  of  finite  width. 

Finally,  I  should  like  to  express  my  special  gratitude  to 
Professor  Betz  for  his  suggestion  of  this  investigation  and  his  active 
support  in  carrying  it  out. 


Translated  by  Mary  L.  Mahler 
National  Advisory  Committee 
for  Aeronautics 
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Figure  1.-    Plate  with  vortex  sheets  starting  from  the  edges  (left),  and 
conEormal  representation  of  the  flow  (right). 


Figure  2.-    Actual  velocity  v    of  a  fluid  particle  and  displacement  velo- 
city  V    corresponding  to  the  similar  magnification. 
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Figure  3.-    Spiral  according  to  the  equation    R  =  H/q)2/3_ 
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Figure  5.-    The  spiral  and  its  image  in  the     ^plane. 
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Figure  6,-    Distribution  density  7(2')    of  the  spiral  in  the    ^-plane. 
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Figure  9.-    First  correction  of  the  form. 
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Figure  12.-    Distribution  density  and  decisive  quantities  after  the 

second  correction. 
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Figure  13.-    Final  form  of  the  spiral  for  the  plate  of  infinite  width. 


i 
0.8 
0.7 
06 

0.5 
04 
Q3 
02 
0.1 

R 

V 

A 

/ 

/T* 

W^Outer  loop 

\ 

tY' 

U      1 

piral   R=-4^ 

^^ 

>::s 

>-a 

y,^/^ 

"    ■■ 

^==^ 

—  - 





y 

* 

1 

<P 

0.577-     77- 


1.57r      ZTT        2-5  TT     ITT        3.57r     ATT 


Figure  14.-    Relation  between  the  radius    R    and  the  angle   cp    for  the 
final  form  and  for  the  initial  spiral. 
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Figure  15.-    Final  distribution  density  in  the  case  of  the  plate  of 

infinite  width. 
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Figure  16.-    Variation  of  the  form  of  the  spiral  with  time  in  the  case  of 
the  plate  of  finite  width  (reduced  to  constant  vortex  magnitude). 
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Figure  17.-    Variation  of  the  form  of  the  spiral  with  time  in  the  case  of 
the  plate  of  finite  width  (actual  scale). 
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Figiire  18.-    Variation  of  the  distribution  density  with  time. 
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Figure  19.-    Growth  of  the  circulation  with  time. 
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